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Abstract. Spiral density wave theory attempts to describe the spiral pattern in spiral 
galaxies in terms of a long-lived wave structure with a constant pattern speed in order 
to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived 
feature. We run an N-body/SPH simulation of a Milky Way-sized barred disk, and find that 
the spiral arms are transient features whose pattern speeds decrease with radius, in such 
a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace 
particle motion around the spiral arms. We show that particles from behind and in front 
of the spiral arm are drawn towards and join the arm. Particles move along the arm in the 
radial direction and we find a clear trend that they migrate toward the outer (inner) radii on 
the trailing (leading) side of the arm. Our simulations demonstrate that at all radii where 
there is a co-rotating spiral arm the particles continue to be accelerated (decelerated) by 
the spiral arm for long periods, which leads to strong migration. 

1 Introduction 

Spiral density wave theory suggested by Q has been the most widely accepted theory of the spiral 
arm for almost 50 years. This theory considers the spiral arm to be the crest of a stellar density wave 
that rotates around the centre of the Galaxy at an angular pattern speed, Qp, that does not vary with 
galacto-centric radius. Spiral density wave theory naturally explains the long-lived nature of the spiral 
arm. However, there are no numerical simulations that reproduce long-lived spiral arms like those 
predicted by spiral density wave theory (e.g. O). Recently, there have been observational studies with 
evidence of a pattern speed that varies with radius, which is diflTerent from the constant pattern speed 
of spiral density wave theory 14111 . We analyse the pattern speed of the spiral arms and the dynamics 
of particles around the arm using an N-body/SPH simulation of a Milky Way- sized barred disk. 



2 Model and Results 

We use our original N-body/SPH chemodynamical simulation code, GCD+ (S). The new version of 
GCD+ includes modeling of powerful supernovae feedback and metal diflTusion of the gas. We have 
set up initially an equilibrium gas and stellar disks embedded in a fixed DM halo potential, and follow 
2 Gyr evolution. The mass and initial scale length of the gas and stellar disks are 10^^ M© and 4x10^^ 
Mq, and 4 and 2.5 kpc respectively. The particle mass for both gas and star is 10^ M©. 

As with the previous studies, we also find that the disk develops transient and recurrent spiral 
arms. We have computed a pattern speed of the spiral arm highlighted in Fig.[Tl and find a declining 
pattern speed with radius (right panel of Fig. [T]), so that pattern and stars co-rotate at all radii, contrary 
to what the density wave theory predicts. Since the arms wind up due to the diflTerential rotation, the 
co-rotating spiral arm is transient (see also |8|). We select a sample of star particles at 5.5-5.75 kpc 
within 20 deg away from the spiral arm and study their motion (Fig.O. Figure highlights the evolution 
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Fig. 1. Face-on view snapshot of stellar {Left) and gas {Middle) disk of the simulated galaxy. The white lines show 
our fitting of a spiral arm. Right: The pattern speed of spiral arm (solid line), and stellar rotation speed (dotted 
line). The pattern speed of the central bar is also shown with dashed line. 




Fig. 2. The normalised density map in the ^ vs. R plane dXt - tQ- AO Myr {Left), t - tQ {Middle) and t - tQ + AO 
Myr {Right). Overplotted dots are selected particles around the spiral arm at ^ = to. Arrows indicate the mean 
velocity field, The size of arrows in the bottom left corner corresponds to 20 km s~^ 



of the extreme positive (red dots) and negative (blue dots) migrators. The 'positive' migrators are the 
particles that migrate towards the outer radii on the trailing side of the spiral arm. They are trapped 
by the potential of the spiral arm, which accelerates them. The co-rotating nature of the spiral arm 
ensures that during migration to outer radii, instead of passing through the spiral arm they remain on 
the trailing side. Therefore, they continue to accelerate until the spiral arm is disrupted. The 'negative' 
migrators are particles that migrate towards the inner radii on the leading side of the spiral arm. They 
are decelerated as they become caught in the potential of the spiral arm, and because of the co-rotating 
nature they continue to decelerate on the leading side of the spiral arm until the spiral arm is disrupted. 
This means that they are accelerated or decelerated for longer, and results in stronger migration along 
the arm. Since angular momentum is exchanged via the co-rotating arm, there is little heating in their 
radial motion |5|. 

We find similar transient co-rotating spiral arms in N-body simulations of non-barred disks and find 
similar stellar motion around the spiral arms lU. We are studying further the spiral arms in our higher 
resolution simulations and endeavour to find observational consequences of the co-rotating spiral arm. 
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